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Abstract 

We have carried out experiments to probe the static and dynamic in- 
terfacial properties of /3-casein monolayers spread at the air-buffer inter- 
face, and analysed these results in the context of models of weak polyam- 
pholytes. Measurements have been made systematically over a wide range 
of ionic strength and pH. In the semi-dilute regime of surface concentration 
a scaling exponent, which can be linked to the degree of chain swelling, 
is found. This shows that at pH close to the isoelectric point, the protein 
is compact. At pH away from the isoelectric pH the protein is extended. 
The transition between compact and extended states is continuous. As a 
function of increasing ionic strength, we observe swelling of the protein 
at the isoelectric pH but contraction of the protein at pH values away 
from it. These behaviours are typical of a those predicted theoretically 
for a weak polyampholyte. Dilational moduli measurements, made as a 
function of surface concentration exhibit maxima that are linked to the 
collapse of hydrophilic regions of the protein into the subphase. Based 
on this data we present a configuration map of the protein configuration 
in the monolayer. These findings are supported by strain (surface pres- 
sure) relaxation measurements and surface quasi-elastic light scattering 
(SQELS) measurements which suggest the existence of loops and tails in 
the subphase at higher surface concentrations. 
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1 Introduction 



A polyampholyte is a polymer which contains both positively and negatively 
charged monomers [Q. A weak polyampholyte is one where the overall charge 
can be adjusted by varying external conditions, usually the pH. This can lead 
to overall neutral polyampholytes, whose net charge is zero, or polyampholytes 
which have a net positive or negative charge. The behaviour of such molecules 
is of interest because amongst their number are found the proteins. An under- 
standing of the way in which polyampholytes behave, in particular the transition 
from an open coil to a compact globule, may give some insight to the under- 
standing of protein folding. 

Polyampholytes are also of interest industrially. Polymers with a net charge, 
polyelectrolytes, are often used as stabilisers, thickeners and in oil recovery. 
However in situations of high salinity, commonly found in oil recovery, their 
thickening effect declines dramatically. Polyampholytes hold out the prospect 
of overcoming this limitation. 

There have been a number of theoretical approaches to analysing polyam- 
pholytes, these include scaling arguments [||, Monte Carlo simulations Q, 
1^, molecular dynamics 0, Flory theory models js) and analogies with 
charged droplets Q. Previous experimental studies on polyampholytes have 
looked at the swelling of gels Q , , or the changes in viscosity of polyam- 
pholyte solutions PI . 

The study of polymers at air-liquid and liquid-liquid interfaces is of techno- 
logical and academic interest , , |Q . This is a result of the importance 
of polymers at interfaces as stabilisers in multiphase systems and as modifiers 
of interfacial properties. More recently attention has turned towards charged 
polymers, which introduce added complexity through electrostatic interactions. 
This is in part through a desire to understand the properties of natural poly- 
mers, such as proteins and polysaccharides, which are usually charged. 

Close to overall neutrality, polyampholytes collapse to a compact configu- 
ration. This is because the chain can rearrange in order to facilitate attractive 
opposite-charge interactions which leads to a compact state. Above a net charge 
eN^/^, where N is the number of monomers and e is the charge density, the chains 
swell. This is because, as the net charge on the polyampholyte increases, it be- 
comes increasingly difficult to rearrange to a configuration that allows opposite 
charge attractions without incurring an energy penalty through the now more 
numerous like charge repulsions. In three dimensional synthetic polyampholyte 
gels the swelling transition is abrupt. As a function of ionic strength, gels with 
a net charge are found to be compact at very low ionic strength, swelling to a 
maximum then collapsing. At very high ionic strength they swell again, this 
arises from the osmotic pressure of a large number of counterions condensed 
onto the polyampholyte. Analogies with charged liquid droplets and Monte 
Carlo simulations suggest that, at overall charges beyond the eN^/^ threshold. 
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polyampholytes will break up into a "string of pearls" , where the number of 
"pearls" is ca. Q/Qc, where Qc is the threshold charge described above and Q 
is the net charge. At high ionic strengths polyampholytes behave like polyelec- 
trolytes through the screening of charge by counterion condensation. 

Measurements of surface pressure (n)-surface concentration (F) curves of in- 
terfacial layers are the two dimensional equivalent of pressure-volume curves and 
in a similar manner can provide information on the inter- and intra-molecular 
interactions controlling the properties of the layer. At low surface concentra- 
tions molecules move independently in the so-called dilute regime. There is a 
surface concentration F = F* beyond which polymer molecules start to over- 
lap, this marks the beginning of the "semi-dilute" regime. The behaviour of 
the n-F isotherm in the semi-dilute regime of a polymer layer can be used to 
measure the Flory-Huggins exponent which is related to the relative 'com- 
pactness' of a polymer molecule. Measurements of interfacial strain relaxation 
can be made which, in common with comparable bulk measurements, can give 
an insight into molecular relaxations. 

In this work in addition to using 'zero frequency' interfacial techniques, we 
have also applied surface quasi-elastic light scattering (SQELS) to probe the 
interfacial properties at high frequency. Langevin [T6| gives an excellent re- 
view of the technique and its applications. Surface quasi-elastic light scattering 
measures the power spectrum of the thermally driven fluctuations of a fluid in- 
terface, and from this information the viscoelastic properties of that interface 
can be determined. The measured thermally driven fluctuations typically have 
wavelengths of the order of 100/im, amplitudes of the order 2A and frequencies 
of the order lOkHz. Here we use SQELS for two particular reasons. Firstly 
SQELS probes the interfacial properties at high frequencies compared to any 
other conventional method, thus it is relevant to rapid events found in processes 
such as emulsification. Secondly we can obtain a dilational viscosity, e', which is 
not otherwise measurable in this frequency regime. This dilational viscosity is 
helpful in validating models for the behaviour of the surface dilational modulus 
with frequency. 

In this work /3-casein was used as a model, weak polyampholyte. /3-casein is 
a milk protein with a random coil structure and a molecular weight of ca. 24kDa. 
It is used extensively in the food industry as an 'emulsifier'. Its emulsifying 
properties are believed to arise from the blocky distribution of hydrophilic and 
hydrophobic residues along its length. We will discuss the amino acid sequence 
in more detail below, but essentially the first 50 residues at the N-terminal end 
of the protein are rather more hydrophilic than the rest of the molecule. 

Neutron reflectometry and ellipsometry experiments |l8| , ||l9| have been 
performed on /3-casein at high surface coverages. They have shown that the pro- 
tein adsorbs to the interface with a thin dense layer right at the surface and 
a thicker, less dense layer beneath it. It has been proposed that this sublayer 
is composed of the hydrophilic N-terminal end of the molecule. This picture is 
supported by proteolytic cleavage experiments |2Q| which show the loss of this 
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terminal region when the adsorbed protein is exposed to a cleavage enzyme. 
Monte Carlo simulations and self-consistent field simulations have been 
carried out using monomer sequences which replicate /3-casein on a coarse grain 
level. These simulations are in agreement with the experimental observations, 
showing a dense surface layer and a less dense sublayer corresponding to hy- 
drophilic regions of the protein. 

The experiments presented here complement these studies because they focus 
on a lower surface concentration, where neutron experiments have not been 
done. Furthermore, we have looked systematically at the behaviour of /3-casein 
on buffers of a wide range of pH and ionic strength. 

Douillard and Agui-Bghin have used scaling arguments to model 
the surface pressure-concentration isotherm of protein layers and interfacialy 
adsorbed multiblock copolymers. These models provide a general framework 
which is not inconsistent with the behaviour of /3-casein, however they do not 
directly account for changes in behaviour with pH and ionic strength. Fainerman 
and Miller have made predictions of the shape of the isotherm at high F 
by considering the thermodynamics of aggregates and their equilibrium in the 
surface layer. 

We have calculatedthe likely charge distribution on the protein from the 



published primary sequence 1 26 and the known dissociation constants for amino 
acids. This calculation shows overall neutrality at around pH=5, in agreement 
with the published isoelectric point An overall positive net charge arises 

below this pH, and an overall negative charge above it. We estimate the net 
charge to be around -10 at pH=9. The total number of charges on the molecule 
is between 30 and 50 throughout the investigated range of pH. We observe 
the well known hydrophilic tail at the N-terminal, and in addition tentatively 
identify two further hydrophilic regions. One of these is more substantial and 
lies in the middle of the chain, the other lies near the C-terminal end. 

There are a number of key parameters important in the study of solutions 
containing charged molecules. The Bjerrum length, 1^, is the distance from a 
charge at which electrostatic and thermal energies are comparable: 

where e is the electron charge, is the relative permittivity of water, eg is the 
permittivity of vacuum, fc^ is the Boltzmann constant and T is the temperature. 
In water at room temperature the Bjerrum length is 7.14A. (Note that defini- 
tions in the literature appear to vary, but this is normally because factors of 47r 
are required in converting between CGS and SI units). The Debye length, k^^, 
gives a measure of the distance at which electrostatic interactions are screened 
out by the presence of ions in solution: 

K-i = {S-kuIb)'^ (2) 
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where, n is the number density of charges. In the NaCl solutions used here, 
r£)=3.035A/I^/^, where I is the molar ionic strength. In this work the Debye 
length varies from 96A at I=0.001M to 3A at 1=1. IM. A final quantity of interest 
is the 'effective' temperature, t. This is given by 

bN 

* " Ib{N+ + N-) 

where N is the number of monomers per chain, and N~ are the numbers of 
positive and negative charges on the chain, b is the statistical segment length 
of the protein. For /3-casein in guanidinium chloride (Gdn-HCl) solution the 
radius of gyration, Rg, measured by neutron scattering is 69 A therefore 
b=11.7A (since Rg=(N/6)^/^b). Clearly we would expect the radius of gyration 
to vary with solution conditions, we simply use these data to show that 1^ and 
b are of the same magnitude. In these experiments ( N++ N~ ) lies in the range 
30-50, therefore the system is always in the 'high effective temperature' limit. 
In order to reach the low temperature limit it is necessary to use a very highly 
charged molecule. 



2 Experimental Methods 

2.1 Materials 

/3-casein (Sigma, C-6905, 90% pure) was used as supplied. Img/ml solutions 
in deionised water were prepared from the dried, powdered protein, stored in 
a refrigerator and used within 5 days. Buffer solutions were made up using 
deionised (Elgastat UHQ, Elga, U.K.) water. Buffers with a range of pH and 
ionic strength were prepared. For pH in the range 5.8 to 8.5 phosphate buffer 
was used, below pH=5.8 citrate buffer was used and above pH=8.5 carbonate 
buffer To control ionic strength NaCl was added, quoted ionic strengths 

also include the contribution of the buffer salts. Buffer pH were measured using 
an electronic meter (ATI Orion, USA) before use. 

2.2 Langmuir trough methods 

Surface pressure vs. area isotherms were measured using a Langmuir trough 
with a filter paper Wilhelmy plate sensor (Nima Technology, U.K.) mounted on 
an active anti- vibration table (Halycion, Germany) , both of which were enclosed 
in a draft proof enclosure. These procedures, which reduce extraneous vibration, 
are necessary in order to carry out the SQELS measurements described below. 
The PTFE trough, area 530cm^ was filled with approximately 500ml of the 
appropriate buffer. Surface pressure - area isotherms of the 'bare' buffer were 
measured before each experiment. The surface was aspirated and the isotherm 



5 



re-measured until it showed no increase in surface pressure on full compression. 
This was to ensure there was no surface contamination prior to the addition 
of the /3-casein. /9-casein was spread on the surface of the buffer by careful 
dropwise addition of the appropriate volume of the Img/ml protein solution. 
The initial volume of added protein solution was chosen such that the surface 
layer was initially in the dilute regime. Typically the volume dispensed was 
ca. 25/il. Lower concentrations of spreading solution were found to lead to 
loss into the subphase on spreading. Repeated spreading was used to explore 
a greater concentration range than possible with a single compression. Since 
hysteresis is observed upon expansion following compression to high pressures 
(n > 10 X 10^^ N/m), no data is presented from compressions of monolayers 
having such a history. These methods has been validated for /3-casein in |]30| . 
The temperature of the trough was held at 22°C by running water from a 
temperature controlled water bath (Haake, Germany) through channels in the 
base of the Langmuir trough. A header tank was used to avoid transmitting 
vibrations to the liquid surface. 
Three types of measurement were made: 

(1) Surface pressure(n)-surface concentration(r) isotherms, the surface pressure 
was measured as the trough barrier was moved at a constant compression rate 
of 40cm^/min. 

(2) Step compression measurements, the surface area, A, was changed by a 
constant amount (AA/A = 0.55) at a barrier speed of 1800cm^/min. The 
subsequent evolution of the surface pressure was measured for 15 minutes. The 
trough area was then further reduced and the next relaxation measured. 

(3) SQELS measurements were made during the relaxation phase of the step 
compression experiments. The magnitude of the pressure relaxation is small 
when compared to the accuracy of the SQELS measurement. Experiments were 
all carried out within 3-4 hours of the protein first being spread on the buffer. 
We do not observe aging effects within this time period. We are aware that over 
a period of many hours and even days, particularly for adsorbed monolayers, 
such aging has been observed. 

2.3 Surface quasi-elastic light scattering 

We use a surface quasi-elastic light scattering (SQELS) apparatus, built in 
house, based on a design proposed by Earnshaw and Hard this is 

illustrated schematically in Fig. |l|. The goal of such an apparatus is to measure 
the power spectrum of light scattered inelastically from the capillary waves at 
the fluid interface as a function of scattering vector, q, measured relative to 
the specular reflection. Photon correlation spectroscopy (PCS) is a convenient 
means by which to measure the small shifts in frequency that this entails. The 
photon correlation is done in heterodyne mode and so it is necessary to provide 
a coherent source of light of the original frequency at the appropriate q value. 
This light is provided using a weak diffraction grating. In order for the het- 
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Figure 1: Schematics of the SQELS apparatus, with inset showing an isometric 
view. Symbols are explained in the text. 



erodyne signal to dominate the correlation function, the ratio of the intensity 
of the inelastically scattered light to the 'reference' light must be adjusted to a 
value of the order of 10^^. 

Turning to Fig. |^: Light, with wavelength of 532nm, is provided by a 150mW 
single mode diode pumped solid state laser (Laser Quantum, U.K.). Polarisa- 
tion and intensity are controlled using the combination of the half wave plate 
(A/2) and prism polariser (P). The beam size, profile and coUimation are con- 
trolled using the spatial filter, S. The grating (G) provides a fan of diffracted 
'reference' beams. The lenses LI (f=150mm) and L2 (f=350mm) perform two 
tasks; they converge the reference beams and the main beam to a single spot 
at the fluid interface and they focus the reference beams and the main beam in 
the front plane of the photomultiplier, situated ca. 2m after the surface. The 
relative intensity of the reference beams is adjusted by inserting a neutral den- 
sity filter (NDF) so that it intercepts the diffracted spots but not the main 
beam. The mirrors M1-M4 direct light from the laser onto the surface and from 
there into the detector. The light is detected using a photomultiplier (PMT) 
and processed using a PC-card based photon correlator (BI9000, Brookhaven 
Instruments, USA), the pulse discriminator used in the PMT is modified to al- 
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low the use of the 'multi-photon' mode originally described by Earnshaw p3| . 
At the detector the laser light appears as a bright central spot with a series of 
focussed reference spots at 2-3mm intervals away from the central spot. Each of 
these spots is composed of the reference beam originating from the diffraction 
grating and inelastically scattered light from the main beam. The reference 
beams are sufficiently weak that inelastic scatter from them can be ignored. 
Each spot corresponds to light being scattered to a different q value, the mirror 
M4 is adjusted in order that the appropriate reference beam falls on the detec- 
tor. The value for the scattering vector, q, and the instrument resolution at 
a particularly reference spot was calibrated by fitting the measured correlation 
function for a pure liquid (water) at that point. Quoted values of the interfacial 
properties are the mean of values derived from fitting 5 correlation functions ac- 
quired consecutively under the same conditions, with each correlation function 
accumulated over two minutes. To reduce stray light, a mirror on the trough 
bottom deflects away from the forward direction any light that had not been 
specularly reflected by the buffer interface. 

SQELS data were acquired from protein decorated air-buffer interfaces main- 
tained in a Langmuir trough as described above. Data for a particular surface 
concentration were acquired with the barrier of the trough stationary, although 
we have found it possible to acquire data from a protein layer undergoing very 
slow compression. 



2.4 Data Analysis Methods 
2.4.1 Dilational Moduli 

The surface pressure, 11, of a protein layer is the difference between the surface 
energy per unit area of the bare buffer, 70, and the surface energy measured 
with the layer in place, 7. Features in the H-T isotherms are seen more clearly 
if the dilational modulus, e^t, is calculated from the isotherm, and this is then 
plotted as a function of T. For an insoluble layer: 

The dilational modulus is the in-plane 'dilational' elasticity of the surface layer. 
For uniaxial stress, as found in the measurements made here, the dilational 
modulus is the sum of the compressional and shear moduli p4| . 



2.4.2 Scaling Exponents 

The XI-F isotherms can be described in the semi-dilute regime using a scaling 
law, with exponent y: 

n (xVy (5) 
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This was introduced by Daoud, Jannik and de Gennes ijlsj, first verified 
for a polymer monolayer by Vilanovc [3^ and has since been applied to a wide 
range of polymer monolayers Douillard j23j and Aguie-Beghin |Q have 

further developed these ideas to apply to multiblock copolymers. The Flory 
scaling exponent, relating the chain radius of gyration, Rg, to the number of 
monomers, N, 

Rg oc N" (6) 

is connected to y by y = 2vl{2v-l) in 2D, and y — 31^/(32^-1) in 3D. The value 
of the exponents are different in 2D and 3D chains, and expected values are 
summarised in Table 1 for different solvent conditions. 
Table 1: 

Values of the two dimensional {v2) and three dimensional (z^3) theoretically cal- 
culated Flory exponents, and values of the corresponding (y2 and y3) exponents 
y (n oc in the semi-dilute regime), for different solvent conditions. 



Conditions 


v2 




y2 


y3 


Extended Chain 


1 


1 


2 


3/2 


Good Solvent 


3/4 


3/5 


3 


9/4 


solvent 


4/7 


1/2 


8 


3 


Poor Solvent 


1/2 


1/3 


oo 


oo 



In 2D y is small ( y=3 ) for good solvents and increases to y=8 as solvent quality 
moves towards the conditions. Here we use the scaling exponent as a measure 
of the overall compactness of the protein rather than trying to link observed 
behaviour to very specific changes. We can measure the scaling exponent from 
the slope of the log- log plot of H-F. However the slope of the e-H plot |2^ also 
recovers the exponent y and this method is preferred because it is not sensitive 
to errors in the amount of spread solution. We find in all conditions that the 
e-H plot is linear in, at least, the surface pressure range — 2 x 10^"^N/m. 

2.4.3 Surface quasi-elastic light scattering 

SQELS data are normally analysed in terms of a model treating the interfacial 
layer as a thin flat elastic sheet at the interface |Q . More recently Buzza et 
al.|Q have proposed a model that explicitly incorporates features of a polymer 
brush into the model of the interfacial layer, introducing bending and coupling 
moduli. Numerical analysis showed that the Buzza model would reduce to a 
thin viscoelastic sheet for the interfacial layer thickness measured for /3-casein. 
We base our presentation of the key results for the analysis of SQELS data 
on the paper by Earnshaw et al |]37| and refer to the paper by Buzza et al 
for some clarification and issues related specifically to polymer monolayers. 
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The dispersion relation D{i-u) for waves at an air-liquid interface, bearing a thin 
viscoelastic layer, is given by 



D{iLj) — [eq^ + iuj?] {q + m) 
where m is 



75^ + iuj?] {q + m) 



pLU 



m 



g2 + i-t-, Reira) > 0, 

7/ 



\iujT] (m — q) ] 



(8) 



T] is the subphase viscosity, p is the subphase density, 7 is the surface tension 
(or transverse modulus) and e is the dilational modulus. 

Solving this equation for D(<j;) = gives us an expression for the wave frequency, 
w, as a function of the scattering vector q. The solutions describe both dilational 
and transverse waves. In a light scattering experiment it is only the transverse 
waves that scatter light and their power spectrum q{}jj) is given by: 



-Im 



iujr]{m + q) 



Diu;) 



(9) 



The behaviour of the dilational waves can be inferred because of their coupling 
to the transverse waves. A fluid-fluid interface can be modelled using a trivial 
modification of (^. In the experiments carried out here a photon correlation 
spectrum is acquired. After accounting for instrumental factors this is simply 
the time Fourier transform, Pg(t) of the power spectrum, P^(a>). The dilational 
modulus can be expanded to take into account viscous effects: 



(10) 



where ep is the dilational modulus and e' is the dilational viscosity. Buzza 
et al have shown that such an expansion is not appropriate for the surface 
tension, and 7' should be set to zero. In this work data is analysed by directly 
fitting the measured correlation function with a correlation function, Pg(t), 
curve calculated from the interfacial properties pTf . The interfacial parameters 
e, e' and 7 are all fitted simultaneously. An example of raw Pq{t) correlation 



data fitted with a modelled function as described above is shown in Fig. 2.4.3 



An alternative approach is to fit the correlation function with a damped co- 
sine which approximates the correlation function calculated using the dispersion 
relation. Then either the fitted frequency and damping or values of the interfa- 
cial properties which are consistent with the values of frequency and damping 
can be quoted |]39| . Since there are three interfacial properties and only two 
parameters it is necessary, in this latter case, to make some assumptions. 



3 Results and Discussion 
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Figure 2: Typical correlation functions obtained by surface quasi-elastic light 
scattering at subphase pH=5.24 and scattering vector q=425cm~^ under dif- 
ferent conditions (a): 'bare' buffer; (b): /3-casein monolayer at concentration 
r = 1 X IQ^'^g/w? and pressure 11 = 6.3 x 10~^A^/m. The solid lines are fits 
with the model described in the text, and insets show the residuals of the fits. 



3.1 Surface pressure (11) - surface concentration (F) isotherms 



Fig. 3.1 shows a selection of surface pressure (H) - surface concentration (F) 
isotherms for /3-casein spread on buffers at a range of pH values, measured us- 
ing a Wilhelmy plate. The accuracy in measuring 11 is 0.02 x 10~'^N/m. For a 
single run the random error in the surface concentration is negligible, comparing 
between isotherms from different runs there is a random error of around 15% 
in surface concentration. However, once this offset in surface concentration is 
corrected for, the isotherms are highly repeatable. 

At low concentrations the pressure is very low, corresponding to a dilute regime 
of isolated proteins at the surface |l^. As the surface concentration increases 
the pressure starts to increase markedly, this is the point where the proteins 
at the surface come into contact with each other and marks the onset of the 
semi-dilute regime F = F* . The behaviour in this regime clearly varies with the 
pH, with the isotherm becoming flatter as the pH, and thus overall charge, is in- 
creased. At the highest surface concentrations the isotherms appear to converge 
to a universal concentrated regime, where there is significant chain overlap. We 
can estimate the radius of a protein molecule on the surface from the upturn 
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Figure 3: Surface pressure (11) - surface concentration (F) isotherms, obtained 
using a Wilhelmy plate, for different pH values, at constant ionic strength 
I=0.01M 



concentration F*. F ranges from 0.25xl0~^g/ni^ for high pH to O.bxlO'^g/rri^ 
at the isoelectric pH. This leads to protein radiuses that increase from 51 A at 
the isoelectric pH to 72A for high pH (due to the uncertainty in the spreading 
procedure the confidence in this result is ±10%). 

Fig. 3J shows a comparison of II-F isotherms obtained using Wilhelmy plate 
methods and surface quasi-elastic light scattering. The error bars shown for 
the SQELS data are the error in the mean for fitting groups of five correla- 
tion functions as described above. At low surface concentrations there is good 
agreement between the static and dynamic measurements. However at higher 
surface concentrations the value of the surface pressure measured using SQELS 
lies around 2x 10~^N/m above the static value. This discrepancy is considerably 
larger than the uncertainty in the fitting of the surface pressure. The onset of 
the deviation always occurs close to the maximum in the dilational modulus, as 
discussed below. 



3.2 Scaling exponent and the semi-dilute regime 



Fig. 3.2 shows the variation of the pressure-concentration scaling exponent, y, 
as a function of pH. Data from experiments with three different ionic strengths 
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Figure 4: Comparison of 'static' H-T isotherms (solid line) to isotherms obtained 
from surface light scattering for different subphase conditions (points): (a): 
pH=5.24, 1=0.01; (b): pH=8.30, 1=0.01; (c): pH=8.34, 1=0.001; (d): pH=7.60, 
1=1.1. The scattering vectors q are in the range 425 < q/cm~^ < 507. Error 
bars are the standard deviation of the mean arising from the average of values 
obtained by fitting 5 correlation functions. 



(I) are included in this figure. It can be seen that the scaling exponent has a 
maximum at the isoelectric point, this indicates that the protein is most com- 
pact at this point. The peak value suggests that /3-casein is close to 9 condition 
and agrees with very recent work on the conformation of polyampholytes |^ . 
As the pH is varied away from the isoelectric point the exponent decreases, 
approximately linearly with the pH. This change in the exponent corresponds 
to a swelling of the polypeptide. The decrease in y is more dramatic at low I, 
indicating that screening charge interactions decreases the swelling effect of a 
net charge on the molecule. The exponent y changes more rapidly on the low 
pH side of the isoelectric point, as expected, since the calculated net charge 
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Figure 5: Scaling exponent y as a function of pH, for buffers with different ionic 
strengths: (A): 0.5 </ < 1.1; (o): 0.008 < / < 0.012; (V): 0.001 < / < 0.003. 
Continuous lines are a guide to the eye for low salt concentration behavior, the 
dashed line for high salt. 



changes more rapidly here than it docs on the high pH side of the isoelectric 
point. Where indicated, error bars are the standard deviation of the mean of 
2-8 repeated compressions, and they indicate the confidence to be expected for 
this kind of measure. 

This behaviour is in contrast to that observed in 3D polyampholyte gels, where 
there is a sharp step in the gel swelling as the overall charge is increased from 
neutrality through the threshold charge of eN^/^. Here we see a more gradual 
change in the degree of swelling. In our system this critical value of net charge 
is anticipated to be ca. 3e. This degree of charging is achieved at a pH very 
close to the isoelectric point, therefore it is possible that the highly collapsed 
state has been missed. Also, the lower dimension of the surface layer compared 
to the bulk gel may play a part in changing the nature of the transition from 
collapsed to swollen states. This dimensional effect has been recently addressed 
in a Monte Carlo simulation of a diblock polyampholyte , where it is argued 
that topological constraints modify the 3D coil-globule transition into a folding 
in 2D. 

Fig. 3.2 shows the variation of the scaling exponent with ionic strength for three 
ranges of pH. 

For pH5-6, the scaling exponent decreases slightly with increasing ionic strength. 
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Figure 6: Scaling exponent y as a function of ionic strength, for buffers with 
different pH: (A): 9 < pH < 10; (o): 7.5 < pH < 8.5; (V): 5 < pH < 6. Lines 
are a guide to the eye. 

corresponding to a small swelling of the protein. This behaviour is consistent 
with that of a polyampholyte, where the screening out of opposite charge inter- 
actions will lead to chain expansion. An alternative explanation for this small 
effect is that increased ionic strength facilitates a small increase in overall charge 
leading to an overall chain expansion due to the increased importance of like 
charge repulsions. 

The exponent y is constant with I for pH7.5-8.5, meaning that the polyam- 
pholyte tendency to contract is balanced by the net charge present on the 
molecule at that pH. 

For pH9-10, the scaling exponent varies considerably with ionic strength. It 
exhibits a minimum at around I=0.003M, i.e. below this value of I it decreases 
with increasing ionic strength and above it it increases. This minimum in the 
scaling exponent corresponds to a max;imum in the chain swelling, and the non 
monotonic behaviour would appear consistent with polyampholyte behaviour at 
low ionic strength followed by polyelectrolyte behaviour at higher ionic strength. 
This conclusion is probably incorrect though, because at the same ionic strength 
as the minimum in the scaling exponent (I=0.003M) we find gross changes in 
the dilational modulus isotherms. Above this ionic strength there are two peaks 
in the dilational modulus and below there is only one. As discussed in detail 
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below, these features probably correspond to gross changes in the surface con- 
figuration of the protein, with parts of the molecule moving into the subphase. 
It is therefore possible that the minimum in the scaling exponent arises from 
this gross structural change rather than a change in behaviour of an essentially 
homogenous series of protein configurations. The distribution of positive and 
negative charges in /^-casein is non-random, it may be that some of the features 
seen are due to this non-randomness. Theoretical work suggests that cor- 
relations in the sequence of positive and negative charges are important when 
the correlation length is of a similar size to the estimated charged blob. For the 
opposite extreme, where the charge sign alternates, the effect is unimportant. 
The value of the scaling exponent at high pH and low salt is in fair agreement 
with the prediction that a polyelectrolyte will behave as an extended chain 
[§. 

Our experimental results can also be compared to previous results on the swelling 
of 3D polampholyte gels as a function of ionic strength, where non-monotonic 
swelling behavior has been observed as a function of ionic strength. English et 
al [lO| found a maximum in gel swelling, Nisato et al ^ see only a minimum in 
swelling as a function of increasing ionic strength. 



3.3 Dilational modulus: e—T isotherms 

The scaling exponents in the previous section probe the proteins when they 
have just started to come into contact with one another. We will now use the 
dilational modulus (e) to probe the behaviour as the protein is further com- 
pressed. Dilational moduli-surface concentration plots for a range of pH and 



ionic strength arc shown in Fig. 3.3. Each figure contains data where one of 
pH and ionic strength is held constant whilst the other is varied, (a) shows 
data for fixed pH = 9.5 and ionic strength varying from 0.001 to 1.1, here two 
peaks are seen in e, one at around F = 0.8 x lO^-^g/m^ and the other at around 
1.5x 10~^g/m^, the magnitude of both of these two peaks increases as a function 
of increasing ionic strength. At the lowest ionic strength the peak at lower F is 
not apparent, (b) shows data for fixed ionic strength=l.l and pH varying from 5 
to 9.5. Once again two peaks are observed in the dilational modulus at F = 0.8 
and 1.5xl0~'^g/m^, although they shift to slightly higher F as pH is reduced. 
In contrast to the data in (a), the magnitude of the peak at F = 0.8 increases 
with reducing pH but that at higher F decreases in magnitude, (c) shows data 
for fixed pH=5.0 and ionic strength varying between 0.001 and 1.1, here only a 
single peak in the dilational modulus is observed which shifts to slightly higher 
F as ionic strength is reduced and shows little change in magnitude. Finally, (d) 
shows data at fixed ionic strength ~ 0.001 and pH varying from 5 to 9.5 once 
again only one peak in the dilational modulus is observed, this shifts to higher 
F with increasing pH and reduces in intensity. 

Note that increasing ionic strength corresponds to reducing the electrostatic 
screening length from 96A to 3A. Increasing pH above 5 leads to increased 
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Figure 7: Dilational modulus Cgt obtained from surface pressure isotherms over 
a range of subphase conditions, (a) and (c) show the effect of varying the ionic 
strength I at fixed pH, (b) and (d) the effect of varying pH at fixed I. 



overall charge on the protein. 

Previously, for /3-casein, only a single peak in the dilational modulus at 
lower r has been commented upon, being attributed either to the collapse of 
the N-terminal end of the protein into the subphase , or to looping of parts 
of the molecule in the subphase |3^] . We propose that the first peak corresponds 
to the tail collapse in the subphase, and that the second peak arises from the 
collapse of a second region of the molecule into the subphase, this may well 
be the loop region identified by examination of the primary sequence discussed 
above. These transitions that depend on concentration, pH and 1 are illustrated 



schematically in Fig. 3.3 



To our knowledge no synthetic system, even a multiblock copolymer, has been 
shown to exhibit two maxima in the dilational modulus as a function of surface 
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Figure 8: Schematic diagram of the /3-casein molecule at the air/buffer inter- 
face. At low concentration F the molecule lies on the surface with a 'pancake' 
configuration. As F increases, first tails then loops of molecule are forced into 
the subphase. 



concentration. 



3.4 Dilational modulus: Surface Transitions 

We have further analysed this data by observing the pressures where maxima 
and minima in the dilational modulus occur, as a function of pH and ionic 
strength (I). By considering a maximum in e as the onset of a conformation 
transition and the successive minimum as the end of the transition, we are in 
the position to present "configuration maps" of the protein monolayer. 



In Fig. |3.4| we show the effect of ionic strength. Shaded regions are for surface 
concentrations where e is decreasing, these are the regions in which transitions 
between configurations occur. The transition corresponding to the first peak in 
e is labeled as t^ (we are proposing that this is the tail protruding in the sub- 
phase), and the transition corresponding to the second peak is ts (this would 
be the looping of some segments). The most obvious effect of salt addition is 
to move all transitions to higher H. This means that transitions become more 



energetically costly the more charge interactions are screened. From Fig. 3.4 



it is now clear how a two peak structure turns into a single peak for very low 



I. At the isoelectric pH, Fig. |3^ (a), one sees that the transition tA does not 
occur at low I. This is consistent with the picture that a neutral polyampholyte 
is compact and tightly folded. Thus there no tail can easily be submerged. On 



the contrary at the other limit of high pH in Fig. 3.4 (c), the tail transition t^ 
appears well defined and at low I it disappears at 11 = 0, meaning that the tail 
is already submerged at F = F* when the polymers first come into contact. 



In Fig. 3.4 we show the effect of pH. The general feature is that as the pH 



increases and the molecule develops a net charge, all transitions occur at lower 
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Figure 9: 'Configuration maps' at fixed pH and varying I, obtained fi'om e.,^-!! 
isotherms. (A) are the location of maxima in est-H plots and (V) correspond to 
minima. Shaded areas, labelled t^ and t_B, guide the eye to regions of decreasing 
e with increasing 11. We tentatively associate t^ as the region where the molecule 
forms tails in the subphase, and t^ as where it additionally forms loops, (a): 
5 < pi? < 6; (b): 7.5 < pH < 8.5; (c): 9 < pi? < 10 



surface pressure. 

In Fig. (a), low I, one sees that the tail transition t^ is not resolved from 
ts at pH5, and emerges from pH around 7.5. It is very thin, meaning that the 
energy difference between conformations with the tail submerged or on the sur- 
face is very small. Also it appears that Ia is moving to 11 = 0, thus confirming 
that at high pH and low I we expect the tail to be spontaneously submerged. 



Moving to Fig. 3.4 (b) and (c) one sees that iA broadens with increasing I, and 
the loop transition t^ shifts towards the upper 11 limit of our study. 
We note that previous studies pO| have shown the /3-casein monolayer to wholly 
collapse into a multilayer at pressures around 22xlO~^N/m, so for 11 > 20 x 
10~^N/m it would not be very significant to explain e behaviour in terms of 
single molecule conformation. 

The concentrations at which the peaks in the dilational modulus occur are 
consistent with the calculated area covered by the protein with, successively, a 
terminal tail and a loop in the subphase. The peaks shift to lower concentrations 
with increased pH, this is consistent with the protein extending as more groups 
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Figure 10: 'Configuration maps' at fixed I and varying pH, obtained from est-II 
isotlierms. (A) are the location of maxima in e^t-II plots and (V) correspond 
to minima. Shaded areas, labelled tA and ts, guide the eye to regions corre- 
sponding to a decreasing e with increasing H. We tentatively associate t^ as 
the region where the molecule forms tails in the subphase, and t^ as where it 
additionally forms loops, (a): 0.001 < / < 0.003; (b): 0.008 < / < 0.012; (c): 
0.5 </ < 1.1 



dissociate. The lower concentration peak reduces in magnitude as pH increases, 
whilst the higher concentration peak increases in magnitude with increasing 
pH. This is not unexpected, the free energy cost of forcing a longer loop into 
the subphase increases with loop length, since the ends of the loop must be 
constrained to the interface. The entropic cost of tail submersion depends less 
on length, as only one end is held at the surface. 



3.5 Dilational modulus: High Frequency 



Fig. 3.5 shows dilational moduli, obtained using SQELS as a function of surface 
concentration, compared with the values obtained from the static H-F isotherm. 
At low surface concentration the static and dynamic values of e are very sim- 
ilar, in common with the surface pressure measurements. At higher surface 
concentrations the dynamic, e, value lies considerably below the static value. 
Once again the divergence of static and dynamic results occurs at or shortly 
before the first major maximum in the dilational modulus. The dynamic value 



20 




Figure 11: Comparison of 'static' e-F isotherms (solid line) to dilational moduli 
obtained from surface light scattering for different subphase conditions (points): 
(a): pH==5.24, 1=0.01; (b): pH=8.30, 1=0.01; (c): pH=8.34, 1=0.001; (d): 
pH=7.60, 1=1.1; 



of the dilational modulus only ever contains one strong maximum, unlike the 
static value which sometimes exhibits a peak at higher surface concentration. 
The significantly reduced value of the dynamic dilational modulus as compared 
to the static value has often been observed with soluble surfactants jl^. In 
these systems the difference arises from diffusion of surfactant out of the surface 
layer during the period of the dilational waves at the surface. This has been 
described theoretically by Lucassen et al and Hennenberg et al Q]. Here 
we do not anticipate a gross diffusional motion of the protein into and out of 
the surface layer at such high surface concentrations. However it is probable 
that the motion of parts of the molecule into and out of the surface layer would 
have a similar signature p5[ . Lucassen et al find, (here we use the notation of 
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Langevin) : 

where ft is the reduced frequency, 




^=\li^.^-\l^^ (12) 



D is the diffusion coefficient of the surface species and c is the bulk concentration 
of adsorbant and thus Tc is the characteristic time for this process of segment 
exchange between surface and subphase. 

3.6 Dilational modulus: SQELS as a function of q at high 
surface concentration 

By varying the scattering vector at which we coUect SQELS data we vary the 
frequency, to, of the capillary wave s we observe. This data can then be tested 
against (pi]). This is done in Fig. 3.6. The data are in the high surface con- 



centration regime where the static dilational modulus is larger than the SQELS 
modulus. A fit to the real part of the dilational modulus was made by vary- 
ing the parameter Tc and the static dilational modulus, e^t- Fitted values are 
e = 9.2 X 10~'^N/m and Tc = 12/xs. The value of est from the compression 
isotherm is 20xlO~'^N/m. If we estimate the value of dc/dF from data pre- 



sented by Graham and Philips |18| the value of D obtained in this way is 
7.5x 10~^^TO^/s. The Lucassen model produces values of the dilational elasticity 
and viscosity of approximately the right magnitude, with the correct dependence 
on frequency, however our data show a value of the complex part (u x e') which 
is larger than the real part (e) and this is cannot be accounted for by The 
discrepancy may be due either to a systematic error in the determination of the 
dilational viscosity or to the fact that the type of behaviour we anticipate is 
simply not modelled well by these expressions. Since the proposed mechanism 
is not simply diffusion into the bulk, the fitted values of Tc might not have a 
real ph ysica l significance. 

Fig. shows values of the dilational viscosity obtained using SQELS. We 
have not measured low frequency values with which to compare this data. How- 
ever we find that the dilational viscosity of /3-casein is very similar to that 
reported in the literature for other polymer monolayers measured by SQELS 
The dilational viscosity exhibits a maximum at the same surface concen- 
tration as the maximum in the dynamic dilational modulus. 

3.7 SQELS as a function of q at low surface concentrations 

In the previous section data from a range of q vectors (and thus frequencies) 
were introduced for layers at high surface concentrations. Here we consider 
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Figure 12: Frequency dependence of the dilational modulus (a) and dilational 
viscosity (b) at pH=7.6, 1=1.1, 11 = 12 x lQ~^N/m. This pressure is in the high 
concentration region, where the high frequency e is consistently lower than the 
static Esf. Solid lines both in (a) and (b) are the result of fitting the data in (a) 
with the model described by eq. (pi]). In (a) the e value plotted at w = is tgt 
from the compression isotherm. 



similar data acquired at low surface concentrations. Fig. 3.7 shows dilational 
modulus and viscosity as a function of q for a protein with surface concentration 
0.8 X 10~^g/m^ on a buffer with pH=8.3 and 1=0. OlM. This corresponds to a 
point close in concentration to the first peak in the dilational modulus, where 
the dynamic values of dilational modulus lie slightly above the static values 
(see Fig. ^^S] (d)). We fit these data using a Maxwell fluid model, i.e. a serial 
combination of a spring and dashpot, characterised by a single relaxation time, 
Tm- The complex dilational modulus is given by: 

/ [i^Tpf \ . / UTp \ 

eo = e^f + eoo— — + M ^°°l~~n ^2 \^^) 

\ I + [uotdY J \ I + [utdY J 

where the first term is the dilational (clastic) modulus and the second part is 
the dilational viscosity, egt is the static dilational modulus, too is the amplitude 
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Figure 13: Dilational viscosities obtained from surface light scattering for 
different subphase conditions: (a): pH=5.24, 1=0.01; (b): pH=8.30, 1=0.01; 
(c): pH=8.34, 1=0.001; (d): pH=7.60, 1=1.1; 



of the relaxation and u is the frequency of measurement. The fit to the data 
is shown by a solid hue, a single relaxation time of 23/is. We believe that this 
relaxation time is related to a tail submersion event during the period of the 
dilational wave. It only appears strongly at high pH and moderately low 1, that 
is, as discussed above (see Fig. ( ^3.4| ) and ( |3.4[ )), the conditions where the energy 
difference between absorbing and desorbing the tail is small. Assuming, as we 
have found, that an expression such as ( [l^ ) holds, then the relaxation time td 
can be obtained from the static dilational modulus and the dynamic dilational 
modulus and viscosity measured at a single frequency: 

TD = — (14) 
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Figure 14: Frequency dependence of the dilational modulus (a) and dilational 
viscosity (b) at pH=8.3, 1=0.01, 11 = 2 x 10~^7V/to. This pressure is close to 
the first conformation transition region, where, for some subphase conditions, 
we observe a high frequency e bigger than the static e^t . Sohd Hues are best fits 
with the model described by eq. (p^). In (a) the e value plotted at w = is e^t 
from the compression isotherm. 



The results of this analysis (data not shown) are consistent with the more com- 
plete frequency measurements. Such relaxation times have been measured using 
SQELS by a number of groups for a number of polymer systems. Monroy et al 
[ pO| found an Arrehnius temperature of the relaxation time for polyvinylactetate 
(PVAc) (M^u = 90000), at the overlap concentration F*. The relaxation time de- 
creased from 80/xs at 1°C to 5/is at 25°C. Monroy also finds relaxation times for 
mixed monolayers of PVAc and poly(4-hydroxystyrene) (P4HS) |Q. The relax- 
ation times are shorter for higher surface coverages and for larger fractions of the 
PVAc. Brown et al |^ find relaxations times of this magnitude which reduce 
as a function of increasing surface concentration for polymethylmethacyrlate - 
poly-4-vinyl pyridine diblock copolymers. Richards et al [|8| find decreasing 
values of the relaxation time at lower concentrations, for PMMA-polycthylcne 
oxide diblock copolymers, but at higher concentrations the relaxation times in- 
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crease with increasing concentration. Mizuno et al |Q also find a Maxwell 
relaxation, and they note that its physical origin can be a potential barrier to 
desorption. 

3.8 Pressure relaxation data 

We have carried out 'step compressions', where the surface layer is compressed 
at a constant rate for a short time and the evolution of the pressure (11) is 
then measured for about 15 minutes, following which the next compression is 
performed. Similar measurements on other systems are reported for example 
by Monroy et al Such measurements have an obvious parallel with step 

strain measurements made on bulk materials. 

We find that a single exponential decay is insufficient to fully describe the data 
(this was also observed in |Q at high H) , but that in every case the spectra can 
be fitted with a sum of two exponentials with well separated time constants: 

n(t) ^ n^e-*/^-^ + nBe-*/^« + Hc (15) 

where 11^ and 11 q are the amplitudes of relaxations with time constants and 
tb- We find that the fast mode timescale is 20 < r^i/s < 50 and the slow mode 
timescale is 300 < tb/ s < 700. lie is a constant pressure offset (the equilibrium 
pressure) and t is the time since the end of the last constant rate compression. 



Fig. 3.8 shows fitted values for these parameters from two sets of step compres- 
sions, measured on buffers at the isoelectric pH with 1=0.1 and at pH=8.3 with 
1=0.01. These buffer conditions arc chosen because they correspond to very 
different isotherms. The ratio between successive areas is always ^A/A — 0.55. 
We are exploring a more complete range of subphase conditions and we refer to 
a future work for a complete presentation of these results Kl| . 



Fig. 3.8 (a) and (c) show the amplitudes Ha and Hb respectively of the fast 
and slow relaxations, as a function of the equilibrium pressure He- One sees 
a general trend for the relaxations to become bigger as the pressure increases. 
At low n negative values of both Ha and Hb are sometimes found. These are 
puzzling, as a negative amplitude reflects the fact that the pressure increases 
after the compression has ended. Although the data presented is too limited to 
support a strong conclusion, we observe that negative amplitudes for the slow 
mode (lis) occur only for buffer conditions that present a "tail submersion" 
transition, and then only around the transition pressure. 

We believe that the fast mode is connected to relaxation of stress in the direc- 
tion parallel to the compression, and that the slow mode is connected to stress 
relaxation between the parallel and normal directions. These stresses arise be- 
cause the fast step compression induces a concentration gradient parallel to the 
compression, and can also bring the monolayer to a transient state where the 
strain in the direction of compression is greater than in the normal direction. 
We further suggest that below the "loop transition" pressure (which is in most 
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Figure 15: Parameters describing the relaxation of pressure after a 'step- 
compression', fitted with the form shown in cq. (p^). (A) correspond to a 
buffer pH=5.60, 1=0.1 and (V) to pH=8.32, 1=0.01 



conditions around 11 = 10 x lO^'^N/m, see Fig. 3.4) the surface pressure relax- 
ations can be understood on the basis of the concavity of the equilibrium II-F 
isotherm. On the contrary, above the loop transition, a fast compression does 
not give time to the monolayer to equilibrate with the subphase, and thus one 
can not refer to the equilibrium isotherm. 

In Fig. 3.S (b) and (d) we show the timescales ta and tb respectively of the fast 
and slow relaxations, as a function of the equilibrium pressure He- 
Noskov ||52j has modelled these relaxations, and his work is relevant to under- 
standing the observed relaxation timescales. He describes the interfacial relax- 
ations as analogous to the Rouse-like relaxation modes of a bulk gel and how 
the presence of loops and tails of the polymer in the subphase leads to some 
degree of entanglement. 
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4 Conclusions 



We have measured the intcrfacial behaviour of /3-casein as a function of pH and 
ionic strength with high frequency hght scattering and conventional methods. 
We find scahng exponents in the "semi-dilute" regime that agree with calcu- 
lations reported in the literature for either polyampholytes or polyelectrolytes, 
depending on the charge present on the protein. In particular, we find that close 
to overall neutrality the protein molecule behaves like a random walk; as the 
pH is moved from the isoelectric pH the protein expands continuously almost 
reaching an extended chain configuration. At high pH and very low salt concen- 
tration, where the protein is overall negative, we observe an effective swelling 
then deswelling of the chain as a function of increasing ionic strength. At higher 
ionic strength the chain acts like a polyelectrolyte. Here only short range inter- 
actions are possible and, since the chain is overall charged, these interactions 
will be predominantly repulsive, weakening as the screening length is further 
reduced thus leading to deswelling of the chain. 

Interfacial relaxation times have also been measured using both surface quasi- 
elastic light scattering and strain relaxation measurements. These suggest that 
at low surface concentrations the proteins act like independent disks, whilst at 
higher surface concentrations the layer forms an entangled network with Rouse- 
like relaxation modes. 

The dilational modulus exhibits one or two peaks as a function of surface concen- 
tration, depending on the buffer conditions. These peaks are associated with the 
collapse of hydrophilic parts of the molecule into the subphase and we present 
a "configuration map" showing the dependance of these surface transitions on 
pH and ionic strength. We suggest, tentatively, that the peak at lower concen- 
trations corresponds to the penetration of a tail from one end of the molecule 
and that at the higher concentration is associated with a loop. This conclusion 
is supported by the discovery of prospective loop and tail regions in the amino 
acid sequence of the protein, whose sizes broadly match the positions at which 
the peaks occur in the dilational modulus. 

The protein /3-casein was chosen in part because as a random coil protein it 
should most resemble a synthetic or model polyampholyte. We believe this is 
the first such study of this type, in applying ideas relating to polyampholytes to 
the understanding of proteins at interfaces. We are currently conducting a sim- 
ilar study into the behavior of a globular protein, for comparison to /3-casein. 
We expect the relevance of the observed behavior of /3-casein to depend on the 
extent to which the globular protein unfolds at the surface. 
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